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REPORT ON NEW CALCULATIONS ON THE MAGNETIC 
POTENTIAL OP THE EARTH 

Adolf Schmidt 
Gotha 


I. In the twelfth ,,r olume of the publication titled "The Archives /j » 
of the German Naval Observatory" , 1 provided some mathematical 
developments on the general theory of the earth's magnetism about 
5 years ago to initiate more precise methods applied up to now 
for studying ti e magnetic state of the earth and also to supply 
some preliminary work for the tasks involved in this undertaking. 

In order to avoid the necessity for research in previous work 
in order to understand the following report, I will present the 
contents as briefly as possible. 

The attempts previously initiated by the work of Gauss for 
representing the magnetic force of the earth with an analytical 
expression have not succeeded in approaching reality, although 
the observation material has gained consistantly in circumference 
and reliability. The differences between the observed a. cal- 
culated values are not essentially smaller than the oldest results 
attained by Gauss himself, even in the newest calculation of the 
magnetic potential of the earth, carried out by Neumayer and 
Petersen. While this gap could be considered a consequence of 
the lack of an empirical basis in older values, however, it must 
now be given real significance and the cause should be looked for 
in insufficient theory, a conclusion first drawn by Mr. Neumayer 
himself from his studies. An improvement in theory is now possible 
and necessary in two directions, without taking into consideration 
the series required for representation. On the one hand, the 
deviation of the earth from a spherical shape must be taken into 
consideration; on the other hand, and above all, the two previous 
assumptions that the magnetic force of the earth exhibits a poten- 
tial and that the origin of this may be found exclusively in 

y. " " ;r ' " 
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the earth's core must be rejected. 


The first and most important task then resulting is to pro- jfk 

vide an analytical representation of the distribution of the mag- 
netic force of the earth on the earth's surface, free of that 
physical hypothesis, a representation which may be calculated 
with a random degree of accuracy, since it depends only on the 
reliability of the observed data and on the expansion of the 
series employed. A sufficient and also most accurate basis is 
provided with the solution of this purely mathematical task for 
all further studies with a final goal of providing a physical 
explanation for the magnetic phenomena of the earth. The question 
can then be considered on whether the entire force effective at 
the earth's surface exhibits a potential or to what degree this 
is not the case, and furthermore, when a potential is found, 
that portion with an origin outside of the earth can be separated 
from the portion originating in the earth's core. 

The flattening of the earth can be taken into consideration 
with little effort and without including the substantial alterations 
in calculation compared to those for a sphere. The series develop- 
ments with a geocentric (not geographical) width or, better still, 
its complement introduced as an argument, , progress as there ac- 
cording to spherical functions? only certain constant factors 
are added with the following significance. A first group of these 
factors depends solely on the flattening of the earth. I have 
calculated this, employing the Bessel number Is a =1 : 299 » 1528 
and publish this in the given location (p.12) under the expression 
p m’ 1T m * ^m the way the s txth order (i.e. for 0 < m 
n < 6) . I compile only the quotients necessary for the present 
purpose irJJspJJ and x^s (Table la and lb) . The numbers in 
these two tables are uncertain by one unit of the sixth decimal 
position, since the values used in this formation, P m * 'H' m etc* 
have been rounded off to 6 digits. A second group of factors 
depends not only on flattening, but also on latitude. These are 
defined by the equations 
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(1) c = Vl -H’^cosv 2 /3 = |/T+7* ^ =s \/ 1 4-* 8 co8t> 8 ; \/ 1 -J-f* 

(t 2 = 0,00671922) 

In this, when a designates the radius of the equator and b the 
polar radius of the earth ellipsoid 


i 2 = (c 8 — 2> 2 ) ; 6* — (2o — 1) ; (a — 1)* 


Furthermore, v is the geocentrical polar distance, connected 
to the geographical, u, by the known equation 


tgt) = v'i+<*tg« 


Table II contains, continuing from 5° to 5°toward u, the appropriate 
values of a , (3), y, as well as the values of a sin v and 
3 sin v, more important for numerical calculations than a and 
3 themselves. In the calculation of v, a greatly abbreviated 
approximation equation was employed; therefore, the indicated 
values deviate from the true values partially by more than 0.5", 
but always less than 1". After I had employed them once for fur- 
ther calculations, a subsequent alteration, requiring a large a- 
mount of correction calculations, was not required because of the 
complete insignificance of those deviations for present purposes. 

I have introduced a deviation from the usual procedure in 

the developments undertaken according to the spherical functions. 

Following the example given by Gauss, the series were ordinarily 

developed according to the functions P n (cos v) cos mA and P 11 

m m 

(cos v ) sin mA with 


m 
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( 2 ) 


P£ (C08 v) = sitt V* [C08 1***""" — C0B 

4 . i!*^3Lfe C0B — , , 

r 2 • 4 • (2 m- — J) (2»* — S) 008 r ‘J 


Meanwhile, the circumstance that these functions are on noticeably 
different orders of magnitude (p|,P^, P^, P^ and therefore also 
P^cos A etc. achieve, for example, the maximum values 1, 0.2 75 1 
0.082, 0 . 023 ) leads to significant disturbances in numerical 
calculations. Above all, the effect of the individual members 
in the series developed according to spherical functions on the 
represented value is not proportional to the coefficients, so 
that it becomes more difficult to rapidly survey the approximate 
course of function and the significance applied to the individual 
members. The coefficients multiplied with spherical functions 
of a higher order, in which m is not completely or almost equal 
to n, are much too large in comparison to their significance with 
the others. Therefore, the coefficients must be applied with 
varying precision since it does not make any sense to employ in- 
dividual members in an aggregate with substantially more precision 
than the others. This uncomfortable and rather raw information^ £ 
since one decimal place already alters the precision of the repres- 
entation rather considerably! has actually hardly ever been employed. 
Instead, as far as I know, all coefficients were calculated to 
the same decimal unit, not justifiable according to the previous 
argument. Thirdly, it should be remembered that the formation 
and resolution of equation systems for calculating the series 
coefficients is facilitated when the factors of this unknown, 
in this case the function values of pJJ, differ as little as 
possible from one another. 


The customary development utilized especially in England 
according to the functions 

(3) mit fi = cos v 


und a* 


ri. 3. 5.. ..(2 >»—!)? 
n\ m ! ’ 


o- = o ' [1-3.fi.. ..(2 h- 1M» 
m (« — m)l 


ffir m > 0 
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-1 


l 






t: 


j m Pn 

<i/r 


sin v* 



leads to the same distortions, only in reverse and more intensive. 

The functions all achieve 1 as maximum value, hut the t!J in- 

' ' tL 

crease very rapidly with increasing n and m in contrast to the P m . 

Por example 


rj = 3P 2 , 21 = 8H, TJs: 10396 PJ 

and the highest values of Tg, Xg, Tg, and Tg are 1.5, 3» 37*1 
and IO 395 . The coefficients of the higher series members are 
therefore much smaller here than they should he according to their 
significance . 

Because of these considerations, X Introduced a certain 
multiple of PJJ in place of it in all developments, chosen 

in such a manner that the square average value of Rjjcos mA 
and Rjjsin mA achieve the same amount on the entire spherical 
surface, specifically 1. (It might appear more practical to 
determine that che average value of R m should always he equal to 
1, where each of the functions Rjjcos mx and Rjj sin mA for m 
greater than 0 would he equal to on the average at a power 

of 2. However, since only practical considerations are to he 
decisive in the introduction of R^, it appeared better to proceed 
in the previously mentioned manner, also advantageous for theory, 
because the members with the factors R^cos mA and Rjjsin mA 
always occur independently in the following numerical represen- 
tations. ) 

Now, when dw designates the surface element of the sphere, /j_ 
in the integration over the entire spherical surface 
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J(p; (cost;) COB miydw = f(P; (cos«;) 6 m W ,A*tf w = _♦* 

( 2 «“HK 


ThiB equation applies generally with respect to tho first integral 
but with respect to the second only when m is not 0. The average 
value of Pjjcos mA and Pjjsin mA is therefore ((2n+l) a . 

since the spherical surface is 4fT. ' It directly follows that r£ 
attains the desired condition when 


(4) R' m (cos v) as f* • P; (cos V) = 1/(2 » -f- 1) a: • P: (cos v) 


The values calculated according to this for the factors r^ are 
compiled in Table III', while Table IV contains the values of the 
functions rJJ from the first to the seventh order 

It should not be overlooked that a complete equality of the 
functions R^exists only with respect to the average amounts. 

The maximum values are different; however, the difference is much 
less than in the case of P**. Under the functions of the first 
seven orders (with the exception of the constant Rq-1) , Rq and 
Rh have the smallest maximum in the amount of 1 . 732 , while the 

X rj 

largest, i.e. that of R^, amounts to 3 *873* The condition that 
the maximum does not decrease with increasing n in contrast to 
the functions P^, but that the increase can be termed favorable, 
because the first, most important members of the development 
are obtained with somewhat more precision than the remaining values 
with the same number of decimal places of the coefficients. 

It is perhaps not superfluous to demonstrate with an example 
the alterations established by the introduction of R^ in place of 
PjJ. For this purpose, I select the development of the magnetic 
potential of the earth for the time 1885.0 according to the cal- 
culations of Neumayer-Petersen, found in the "introductory remarks" 
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in the fourth section of the Berghaus Physical Atlas (p. 19)* 
According to the remarks on page 5 of my report mentioned at the 
beginning, I provide all numbers with the opposite sine* moreover, 
I introduce the quantity O.l^cnf" s’" 1 . With these alterations, 

I take the following coefficients from the indicated publication 
for the series for (V:R) according to the function pjj# 


m 

1 

2 

8 

4 

9*o 

—31572,0 

— 790,6 

2436,3 

3439,5 

9* 

— 2481,4 

4979,8 

—3956,0 

3059,7 

K 

6025,8 

—1299,9 

— 738,3 

1187,7 

9s 


566,7 

2785,7 

1975,4 

h; 


1260,4 

44,3 

— 714,7 

9% 



327,0 

— 684,2 

K 



649,2 

— 512,1 

9*4 




84,9 





— 96,8 


The coefficients of the development according to the functions 

R?, on the other hand, are 
m 


n: 

1 

2 

8 

4 

9* 

—18228,1 

— 235,7 

368,3 

262,1 

9*i 

— 1432,6 

1285,8 

— 488,3 

184,3 

»? 

3479,1 

— 335,7 

— 91,1 

71,6' 

9*2 


292,6 

543,7 

168,2 

K 


650,9 

8,6 

— 60,9 

9*2 



156,3 

— 109,0 

K 



262,6 

— 81,6 

9*4 ' 




38,3 

K 




— 43,7 
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After these preliminary explanations provided to avoid a 
disruptive search for references for the following presentation! 

I will present the course of calculations after my earlier work 
in the shortest possible form. In this work, I still assumed 
the development according to the functions meanwhile, I 

simply transmit the designation* introduced therr to the corre- 
sponding developments according to the functions R 5 *. The modifi- 
cations therefore necessary in some equations, limited to the 
entrance of constant factors dependent on rjj, are easily seem 
Therefore, I believe that I can introduce them without special 
emphasis, just as some less important alterations in the design- 
ation. 

The observed values of the magnetic elements of the earth j_ 
at the most numerous possible points of the earth's surface form 
the empirical basis for the entire calculation. First, the com- 
ponents of force X (horizontal to the north) , Y (horizontal to 
the east) and Z (vertically downward) are to be calculated from 
them, replaced by the quantities aX 0Y y Z, termed X, Y, Z, de- 
viating only slightly from the first values with respect to the 
spherical shape of the earth. 

The first task is now to find an analytical representation 
for the distribution of the values of "X, T, Z over the entire 
surface of the earth. This task ran be solved, considered theo- 
retically, with randomly precise approximation of the true condition 
by presenting the quantities Z sin v, Y sin v and Z by rows, 
progressing according to spherical functions oi the arguments 
v and A ( the geographical length). (In the case of I and 1, this 
form of representation is not possible, because it becomes variable 
at the poles. In contrast, however, other expressions, for example, 
(X cosA +Y cos v sinA ) and (X sinA ~Y cos v cos;\ ) , occur in 
place of X sin v and Y sin v.) 


ORIGINAL WW.Z ft 
OF POOR QUALITY 

From this series found for X sin v and Y sin v, the functions 
are now derived 


u- Jxdv t 

o 


X 

W = i)i {v) — J F sin t> 

o 


where ip (v) designates the portion of U independent of X . When 
both are identical , they represent the potential V of magnetism 
for the earth on the earth’s surface with the exception of a con- 
stant factor? specifically, V is the common value of bU and bW , 
where b signifies the polar radius of the earth. 

When U and W differ, this difference provides the proof 
that the magnetic force in the earth’s surface does not exhibit 
any potential, leading to the conclusion of the existence of elec- 
trical . currents, passing perpendicularly through this surface. 

The density, i.e,, the intensity of these currents in relation 
to the surface unit, results clearly from the difference (W-Tj). 

The manner in which these currents are enclosed inside and outside 
of the earth's surface, however, remains completely undetermined 
and therefore the portion involved in the values of U and W, 
therefore also of X and Y and those of Z. 

If it is now required that this portion is a minimum, /1 0 : 

especially that Z is never altered, a large portion remains to 
which a potential is attributed. This is again designated by V. 

Finally, when the series development of V is combined with 
that of the vertical components "Z , the potential stemming from the 
external and internal agents can be separated. 

The final result is therefore formed by three functions 
V i( V a and i sin v (~ i sin v) , precisely expressing the con- 1 

dition of the magnetic field on the earth's surface represented 
in the series for X sin v, Y sin v and IS, therefore also randomly 
approaching the true values through completion of the observation ' 
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data and continuation of the series development, is the poten- 
tial of magnetic masses or enclosed currents in the earth’s core, 
V a is that of fauch agents in the external space , finally i 1 b the 
intensity of the current passing through the earth's surface 
perpendicularly, to he comprehended inside and outside as enclosed, 
so that its effect in this surface itself is a minimum. 


The course of the calculations described' in the previous 
section is represented by the following equations. 


X sin V ss 2 n n m ( B* m cos mX 4 - <7; tin tnX) 
( 5 ) f sin v = -i" JR* (jD; coa mX 4 - ein mX ) 

Z = 2. (j* cos t»X -f- K sin mX) 


The symbol - n? stands here and in the following for ZoZo 
The development of 1 is only subject to the condition that j° 
must be equal to 0/ in contrast, those of 1 sin v and T sin v 
are connected by a series of conditional equations with the final 
expression that the horizontal component of the magnetic force 
of the earth in addition to all the differential quotients at 
both poles (i.e. for v-0 and v-l80°) is finite and has a clearly 
determined direction, ,i#e. that the poles are not variable points 
for them. These conditional equations are, when near the north- 
pole (v=0) is equal to sin v ra with the exception of infinitely 
small quantities of an higher order, therefore near the southpole 
(v= 180 °) equal to (~1 ) n ~ m ajj sin v m and when 


= r: + fz+'+..* 

4 - t'’:*' + . . 



is intro due edt 
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J?- js: = o 


0: + &a* m & m ~ 0 

r 1 ai c: + ^ «: z>: « o 


When m=0, that sum of the two disappears from which the 
left side of that equation is composed, since She function 

values ajj are determined by the equation 


( 7 ) •: = 2 ~ =$£$ k 

Especially 

a” = rX, <C +1 = r* +1 , aX +i = rX +t , a"+* — 2m ~^ , 2 ,-*<+• 

’ " ’ 2wi*f-8 r ’ " 2wr^5 *■ 


(The corresponding data on the values of PJJ- for v«0, mentioned 
in the above quoted article p. 21, and the equation (14) drawn 
from these are incorrect with the exception of m=0.) 

(The method developed here for the analytical representation 
of the components of the earth's magnetism can of course be applied 
to all cases in which a constant function (a vector function) 
is present, defined clearly everywhere in size and direction, 
while the usual, known development is employed for a variable 
scaled quantity. When the direction of the vector falls into the 
ellipsoid surface everywhere, the developments given in relation 
to X and Y are sufficient; otherwise, a simple representation 
of the vertical components according to spherical functions is 
added as here, generally without the condition Oq-0.) 
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The immediate result of the development found for Y sin v 
is 


W — V' (fO *— P3+ ^o-®o + tyRl -f m.)^ -j- -2'JR* (j~ jEJ cob mk — D^Bintuk) 
( 6 ) ~ V(v) — sin t? (j) (t/) l 4- Wj ( m > 0) 


The derivation of U from X sin v is not as simple. When 


J sin tT" 1 dv = TJ n 

o 


is set and a series of new coefficients ” m mX sin »»*) , 

GjJ, hJJ is introduced through two equation systems provided in the 
following, the result is 


111 


( 9 ) 


U _ (n, IJ, -j- n t I7 t -j- ) -j- -ijR" (£?* cob mi -f- E* m sin mk) r= f(v, k) -f- £7 0 


In this case especially 


¥'(») = G&+GJ.RJ + <??£? + .... 


In order to obtain the above-mentioned equation systems in 
a form most suitable for the numerical calculation, I introduced 
the following designations: 
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( 10 ) 


*i = l, 


3P — 5P-2 («»+jl)(w + P»-lL 
" ~ " m+JJ-1 


( m +p— 1)^^+^* «i* 


In order to avoid applying the equations twice (for the co- 
efficients of the cosine and for that of the sine of mX ), I also 
write 


B n m COB mX 4- C* m Bin mX — Al, O n m cos wX + H n m Bin viX = FZ 


Now the equations are expressed as follows: 


( 11 ) 


7 i m = </° AZ 
J'i jFH +1 — 
j/j, i^+» = • 


^ +, 4-/^: +4 -b 
aU : +4 + 


( 12 ) 


r l m FZ = /4 <K +X 




It can be seen from these equations that Y sin v must be 
developed just as far, but X sin v must be developed up to the 
members of the next higher order in order to obtain the development 
of U and W up to functions of a defined equal order (n). 
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After U and W have been calculated, on the one hand, the 
potential results 


(13) 



(Po + ^o) = i 2& m 0; COB mX -f hi Bin mX) 


on the other hand, the current density of the flow passing through 
the earth's surface vertically 

(14) i = 1- 3'(W- U) 

4fta(ibBiBv dvdX 


The currents directed into the earth are designated positive 
in this case. 

Finally, the coefficients of the potential of inner forces 
result 


(15) 


F, ~ i ^Jtl (cl cob mX ~j- si Bin tnX) 


and those originating from external forces 

(16) F a — (yl cob niX -f- sin mX) 


through the equations 


0F P0 ^ QUALITY 




( 17 ) 


c - - o: — Jzj: *: = c hi — <% % 

r- ~ & — c a* = h n m — si 


in which 


(18) 


J " = + (»+i)-!-) 


und 



d; 


are constants dependent on the flattening of the earth. In Tables 
Ic and Id the values derived from the numbers of la and lb are 
compiled (the coefficients c,s have been mistakenly exchanged 
with y a on page 23 of the above-mentioned work in the equations 
presented above). 

The final result of the entire calculations are the three 
series progressing according spherical functions for V.:b, V :b 
and a|3 bi, with the later remaining behind the two others by 
one degree in the order of the highest members of the development, 
as can be easily seen. 

From these two series with coefficients completely independent 
of one another the series for X sin v, Y sin v and Z can be 
extracted in reverse when considering the regulation established 
for the derivation of V:b from U and W. This requires no further 
explanation. To the above compilation of equations, I will add /1 4 
some developments which I noted in the work mentioned above without 
deriving them, so that I will deal with this briefly here. As 
is known, it is possible to assume a clearly defined distribution 
of free magnetism on an internal surface infinitely adjacent to 
the earth's surface, corresponding within the earth's surface itself 
and in the entire external space to the same magnetic effects, 
as are exerted by the actually present agents recognized as cause 
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of V^j furthermore, it is possible to replace the external forces 
with the potential V 0 in the earth's surface by defined magnetiz- 
ation of an external surface infinitely close to it, insofar as 
its effects are noteworthy in the innerspace and at the surface. 

More important because it is probably of greater significance 
for the physical explanation of the earth's magnetism is that the 
same success can be achieved through an arrangement of electrical 
currents, as is known, filling an infinite adjacent surface parallel 
to the earth's surface (in reality, instead of the surface a 
layer of finite thickness must be assumed as in the first case; 
however, it is easy to see that when the thickness and distance 
from the earth's surface are only small in comparison to the 
radius of the earth, the condensation of the entire current caused 
by simple vertical shifting with suitable reduction in density 
alters the total effect in an area only by small amounts of higher 
order.) The distribution of these currents and the magnetic 
occupation will now be calculated. 

The potential of the free magnetism spread in an area (of 
never infinite large curvature) amounts to the same at every point 
on both sides of the surface. More precisely, when V £ designates 
the potential at a point with a distance from the surface e cal- 
culated positively from a defined direction, V and V have 

+ 8 —S 

the same limit value within one and the same normal with infinitely 
decreasing e. This is represented in a frequently applied sym- 
bolic manner 


r+o ~~ 


F-o 


The alteration in potential or, practically the same, the 
component of force perpendicular to the surface, in contrast, 
is the same in absolute amount on both sides, insofar as it stems 
from the infinite adjacent portions of the magnetic field, but 
in opposing directions. The difference of the values on both sides 
depends on the surface density p at the appropriate position. 
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Concerning the effect of the remaining field situated at a finite 
distance, it supplies the amount of 0 to every difference as 
constant variable* therefore, this remains unaltered and is 
dependent upon p # Written symbolically this is 


(1Z) -(H) _ , 

\StiUo \dn/-o 4:7 P 


or in the application to be considered here for the magnetic force 
in the earth's surface 


Z + o - = — 4»c» 


In this case it applies that the direction of the normal value 
from the inside to the outside, i.e. from below to above, is pos- 
itive, the reverse situation from Z. 

On the other hand, when the simple, whole surface is subjected 
to electric current, the magnetic force caused by this is constantly 
variable when passing through the surface* especially, it is also 


pi) =(lL) 

V Sn /+c v 3n /— o 


or 


J +o 
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On the other hand, the potential at this point is variable. 
The difference of the values at both sides in infinitely adjacent 
points of the same normal value is proportional to the total in- 
tensity S of that current passing in the surface, including the 
observed point* 




— V _ 0 


Those currents passing in a positive direction (counter-clock- 
wise) considered from the side of the positive normal value are 
set as positive in the calculation of S and the others are set as 
negative . 

The total intensity S is now apparently the same for all points 
of a current line and it is altered in the transition from such 
a line to another one by the amount of the current flowing between 
both. It is therefore 

S-constant 

the general equation of the current lines, and the total current 
present between S=S.j and S=S£ has the intensity Sg-S-j . The current 

density is inversely proportional to the alternating distances 
of the lines S=constant, and the direction of current is such 
increasing values of S on the left when viewed from the side of 
the positive normal values. 

In the following equations, the developments serving for 
the definition of the magnetic condition of the earth’s surface 
are introduced, insofar as this exhibits a potential. This is 
according to equation ( 15 ) and (16) 

V, - bSR* (C cos ml -f s* sin ml) V 0 = b cos ml -f d* m sin ml) 
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or more briefly 
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V t — bSJI 


r.^bsji 


can easily be continued in the external space, since it 
stems from internal forces, V into the earth's core for corresponding 

cl 

reasons. However, V. here and V„ there are first not determined. 
Accordingly, I set 


(19) 


(F,)+, = b2Jlfy~ m 

(v,u = 


(V.)-. = 


In order to simplify the equations, I will drop the continuously 
repeated indices m and m everywhere in the following (also in 
the newly introduced quantities M, N, tt , p, x> <1» <$ , e)» and 

the result will be no more disturbing than the circumstance that 
the functions A and B m introduced here temporarily are in agreement 
in the expression with the coefficients occuring at another point 
for X sin v. 

In order to derive Z from V, I must refer to equations (18) 
and (21) c£ my earlier work, according to which 




the value for Z results 
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z = TJ = -H- (n+1) 7 Af+ "7* ) 

Therefore, in the present case 


(%<= -7-s<«+i)-j-J r 


(Z.).,= 


When v + o =V -0’ order represent the observed forces 
as the effects of a magnetic surface field on the earth, it fol- 
lows from (19) (for all indices n, m) 

J=K A~B 

In retrospect to (18) and (20) at 


4:71 = (Z ( ) +Q _ (£,)_<, 



~ 4 71 ( J . = (3,) +0 — ( ZJ _ 0 



On the other hand, when it is necessary to attribute the 
observed effects to currents identical to Z +Q =Z the equations 
are derived from equations (20) 


.A** ) 3f 


or poor 



— (n+l)"- J == — n^A 


and the result from (19) and (18) with a slight conversions 


4 ”%~(Kko-(K)„ 0 

e 


4 * s '*° (r.) + o~(Ku 



In the equations gained in this manner, expressed as follows 
in a detailed representation 


(21) fc-j ~zi m K(r:™ m x+o:'in ml) 

(22) g, = ± cob *»*+£ ein mty, S,=-j 


the complete solution for the two tasks presented at the beginning 
of these developments is found. 


XX. From the beginning of my occupation with the studies 
presented in the previous section, it has been my understandable 
wish to apply the theoretical developments to the actual state 
of the earth’s magnetic force. I was strengthened in this desire 
by the results of an initial experiment undertaken in 1886, be- 
cause the observations did not appear to agree with the hypothesis 
of a potential for the only portion of the force under consideration, 
that is the horizontal portion, but a sudden decision on the chosen 
basis was impossible because of the relative slight difference 
in the contradictions. I had based the calculations on the charts 
of declination and horizontal intensity for 1880.0 published by 
the German Naval Observatory and already employed previously 
by Quintus Xcilius for a determination of potential, charts with 
a scale, not permitting an expansion in the removal of the functional 
values beyond 0.001 cm" : %^ s~^at H and 0.1 at <f . When I found 
time several years later to reassume the task, it had to be my 
desire to use more precise values for the final calculation. 

The difficult and time consuming work of obtaining these with suf- 
ficient completeness became unnecessary because of a very lucky 
circumstance, also with respect to the factual basis. Director 
Neumayer wanted to subject the potential calculations carried out 
by him together with H. Petersen to a controlled repetition and 
he requested that I replace his honored associate after the death 
of Mr. Petersen. I did not delay a moment in accepting this re- 
quest, providing me with material"^, unique in completeness and 
critical research — processed by a researcher in the field 
for decades in scheduled work with all means of personal influence 
and collected at a central position such as the German Naval Ob- 
servatory and prepared with a great deal of experience. 

n r — — 

'Precise information is supplied on this by the text "Atlas 
des Erdmagnetismus" (Department IV of the Physical Atlas by Berg- 
haus) o Furthermore, Dr. Neumayer *'s lecture at the VIII German 
Conference in Berlin (1889) supplies a good survey, "on the pres- 
ently available material for research on earth and world magnetism." 
(Verh d. VIII. d, Geogr.-T.zu Berlin, p. 33ff») 


Mr. Neumayer expressed agreement in my carrying out the new cal- /19 
culations according to the completed theory, providing sufficient 
indirect check of the earlier calculations as a side result. 

Mr. Neumayer did not limit himself in supplying me the basis 
for my study i I also thank him and X am glad to he able to 

present this thanks publically -- in other respects for ener- 
getic support for my work. His lively interest in its progress 
led to an almost uninterrupted scientific correspondence, providing 
me with good suggestions and the joy of critically founded agree- 
ment and encouragement. His repeated support also made it possible 
for me to have a considerable portion of the numerical calculations 
carried out by an auxiliary calculator. My father, who had already 
performed valuable and reliable help on earlier occasions in this 
manner, assumed large sections of the complicated calculations, 
probably not mastered by most calculators* and X would also like 
to thank him at this point, so that I hardly had to carry out 
half of the calculations in the end, The conclusion of.. -the 
work would have been even more delayed than was already the case 
without his assistance. 

After I had received the basis for calculation in the summer 
of 1892, to be reported further (p. 21, 22), I first carried out 
a preliminary development, limited as were all previously publi- 
cized to the members of the four first orders. The results of 
this work, concluded in the spring of I893, have already been re- 
ported briefly in various publications. At the fall meeting of 
the "Association of Friends of Astronomy and Cosmic physics" 
in Muenster (1893) » Mr. Neumayer provided some information on 
this, after I had briefly done this for the physical department 
of the natural science meeting in Nurnberg at their wish. He had 
undertaken a more detailed report earlier in an article sent by 
him to the meteorological congress in Chicago. This article 
has been lost in the meantime in an unexplainable manner. 
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The following pages contain the report on the expanded rep- 
resentation proceeding to the members of the sixth order, carried 
out since then* Because of the condition that this was an initial 
attempt to find a now way, the main emphasis in the discussion 
of results has been placed on the question of the tasks resulting 
for further research and the conditions to be fulfilled in order 
to promote this aim. When certain practical demands are ascer- 
tained as unavoidable in this case, already emphasized for other 
reasons, it may be permissible to express the hope he^e that this 
result will contribute to the acceleration of fulfilling these 
demands. 

In the numerical calculations, which I will now discuss, 
it was above all the aspect of facilitating further repetitions 
with altered empirical basis the determining factor. For this 
purpose, all operations were carried out generally, insofar as 
this was possible, and the observation data employed as late as 
possible, (this produces the further advantage that the general 
solutions derived in this manner can also be applied for treating 
other geophysical problems).. Furthermore, for the same purpose 
the calculation of the values of the components of force following 
from the analytical derivations were examined for a large number 
of points, work which only has been completed to a small portion. 
After this is concluded, it will be sufficient on the basis of 
new material in repetitions and improvements in the work to intro- 
duce the differences of the observed values in comparison to those 
already described and calculated analytically into the calculations. 

E. Sobering has already described in detail the best method for 
l ) 

this purpose. ' 

My initial intention of publishing all results and the complete 
course of calculation in order to facilitate a detailed check 
or possible alterations had to be given up for the time being. 

Tf " 1 1 

Compare K. Schering, "The progress in our knowledge of magnet- 
ism in the earth," Geogr. Jahrbuch XV, 1891, p. 143 ff. 
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I had to limit myself to including first the basis for calculation 
and the most important results? meanwhile, I hope to find the 
possibility for a more detailed publication at a later time. 

The empirical basis for the derivations in the following 
is the same, as was already noted, as that for the calculation 
of potential by Neumayer and Petersen. This consists of the 
trigonometrical series representations of the components of force 
X, Y, Z for the 25 equidistant parallel circles of 0°, 5°> 10° ... 
60 °of the northern and southern latitude, continued up to the mem- 
bers of the fourth order. The coefficients of each of these 75 
series are based on the 7 2 functional values, occuring according 
to the original charts of Dr. Neumayer in the points of intersec- 
tion of the corresponding parallel with the merideans of 0°, 

5°, 10° . . . , . 355° eastern longitude of Greenwich. Therefore, the 
analytical representation of each component is based only on the 
values in 1800 points, considered as observed, distributed rather 
evenly over the earth's surface, although excluding the polar 
areas. In order to provide a view of these values, I compiled 
some of them in table V. 

The series are applied according to the designations of 
Gauss in the form 


*0 + *1 cos;, 4 - Ki sin 4- -f shU/t 

r = l 6 -j- Zj COS X L x Bin * -f 4- z , 4 Bin 4 * 

Z = m 0 4- vtj COB 4- Jf, sin i 4- ...... 4. sin 4 X 


Their coefficients are compiled in the tables VI a, b, c. The 
unit employed here, as in all further data is 


0 , 1 ® cm-i s~ 1 


i.e. the 10000th portion of the Gaussian unit. 
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It is now necessary to develop these coefficients or rather 
(p. 9) the products 

cd* ein t>, aK* sin v\ ftL' m n in t>; yrn* m , yMl 


according to spherical functions. The distribution of the given 
values on randomly selected parallel circles does not make it 
impossible to apply the Neumann method 1 ^ to the development accor- 
ding to these functions, but it is disadvantageous, since the series 
terminates at a large distance from the poles, In order to avoid 
this, and to apply that method in its simplist form, I would have 
had to employ the cartographic representation of the elements 
and determine their values at a number of new latitudes chosen 
correspondingly, In fact, I consider it necessary to carry out 
this study subsequently. In the present case, an important ad- 
vantage would have been lost in this manner. The possibility 
is gained of observing the effect of the altered manner of treat- 
ment very precisely, because the calculations carried out by 
Neumayer and Petersen according to the simplified theory are /22 

based on precisely the same empirical numerical values as those 
undertaken here according to the completed theory. Because of 
the considerable differences between observation and calculation, 
on the other hand, the comparison of results obtained in both 
cases would have suffered in the application of different parallel 
circles in one and another case in a manner not to be overlooked 
when the cartographical representation used as a basis had remained 
the same, as in this case. 

"*■) Astronomische Nachrichten, . Vol. » 15, , p. 313 and Mathematische 
Annalen, Vol. 14, p. 5^7 (reprint). I have given a short compil- 
ation of equations (without explanation) in my first article 
(Aus dem Archiv der Deutschen Seewart, 1889, no. 3). Seeliger 
provided a more detailed discription with tables in the session 
reports of the mathematical-physical class of the academy of 
Sciences in Munich, 1890, p. 499. 
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Therefore, I have employed the method of smallest squares 
for deriving the development in the customary manner, f ^ is 
any of the coefficients 

a f hm,< sintv, a ( K m i sin v„ in«„ /9, , sin ©„ y t n y i M m>i 

where the index i passing through the values 1, 2... 25 differen- 
tiates the values belonging to the different parallel circles 

w, = 30°, u 2 = 35°, = 150° 

and F m designates the corresponding value under the coefficients 

C . D . E . j. k . The system of error equations is then 
nr m m m' °nr m ^ 

expressed as 


n=m+r 

fm,< = Zj K • K (C08 V t ) i = 1, 2 , . . 25 

US zm 


The resolution of this system also depends on the weight 
determination. It could be possible to introduce the equations 
refering to relatively unknown areas, e.g. to higher southern 
latitudes, with less weight then the others. It is, however, 
easily seen that this entails the danger of buying the gain in 
reduction of random error with a systematic distortion of the re- 
sults in this manner. Moreover, the question must be considered 
from another aspect. The quantities for which the series develop- 
ment is sought are, when the flattening of the earth is neglected 
for reasons of simplicity, X sin u, Y sin u and Z, while the 
observed values, are those of X, Y, Z. Compensation would then 
be undertaken by the most precise representation of the latter, 
while this is carried out for the former in the assumption of the 
same weight for the error equations. Accordingly, the weight 
should only be introduced as constant in the development of Z , 
while it should be set inversely proportional to the square of the 
sine of the width in that for X sin u and Y sin u. On the other 


hand, it must he considered that the areas around the equator 
actually do not have such densely arranged, observed values in 
the calculation as the areas closer to the poles because of the 
decreasing length of the parallel circles with increasing latitude, 
of course, disregarding the exceptional polar areas. In order to 
compensate for this differing density, in other words, to permit 
an effect on the result on the basis of the individual zones 
corresponding to their area, it would apparently lie necessary 
to multiply the weight of each equation with sin u in all three 
components. The combiiiation of these different aspects would 
therefore lead to the setting of the weight with (It sin u) in 
the development of X sin u and Y sin u, but in the case of Z 
to select the reverse sin u for this purpose. 

When it is now considered, however, that in the relatively 
large number of excess equations the selection of weight has only 
a slight effect on the result, insofar as this stems from random 
errors, since in this case the ratio of the extreme values for 
weight is only that of Is 2, it appears justified to decide the 
question according to practical considerations. These, however, 
support giving all equations the same weight. The establishment 
of normal equations becomes somewhat easier then and, above all, 
the resolution is the same for all three components, while two 
systems of normal equations — one for X sin u and Y sin u, 
the other for Z — otherwise have to be prepared and solved, 
Finally, a justification based on fact can be supplied for this 
procedure. The present task does not finally deal with the compen- 
sation of the observed state considered erroneous according to at 
theoretically d.erived, enclosed equation, but rather the analytical 
representation of the state as it has been observed in an arbi- 
trarily cut-off derivation, therefore only approximate. Therefore, 
as long as the remaining differences are not reduced below the 
amount of suspected uncertainty of observations, they are to be 
decreased by continuation of the derivation. The task then at 
least becomes approximately defined and the determination of ^weight 
is no longer important. This situation becomes" most clear when 
the observations are considered as distributed in sufficient density 
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over the entire surface of the earth. In this case, where the 
individual coefficients are represented independently of one another 
by the known integrals, the task of compensation is relegated to /2 4 

the background. It is apparently appropriate for the procedure 
to be applied then, when Z is treated in the same manner as X sin u 
and Y sin u or to return to the more precise representation 
yZ, as well as ax sin v and 3Y sin v in accordance with these 
considerations. I have introduced all equations with the same weight 
in the formation of the normal equations. 

Concerning the limitation of the derivation, the final de- 
cision on this mater naturally depends on the permissible deviation 
between observed and calculated values. For the present, it ap- 
pearedpractical to undertake an attempt with a derivation extended 
to the members of the sixth order, since the members of the first 
four orders are not sufficient without a doubt according to previous 
experience, on the one hand, and, on the other hand, the extent 
of the calculation, rapidly increasing with growing number of members, 
forces a limitation. (Proceeding only one order further would 
have been half a step, because the members of even and uneven 
orders are found independently of one another as a result of the 
symmetrical distribution of the 25 parallel circles employed.) 

Formally, it must be required of the limitation for a deri- 
vation that it is not altered by a mere transformation of coordin- 
ates. It would follow from this, as the derived expression for the 
for the spherical function P n (cos v : cob sin v l sin v 2 cos (A, — ?->)) shows: 

that the series must conclude with those members having agreement 
for the indices m and n. Since in the method employed here as 
previously, first leading to the derivation of individual parallel 
circles, obtained independent of one another from the various 
multiples of the geographical length of the dependent portions 
of the function to be derived, there is no objection to defining 
the maximum value of m without consideration of the value of n 
with the exception of the understood condition that in each indi- 
vidual member m may not be larger than n. 
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I have expanded the derivation up to the values m-4, n-6 
(in the case of ax sin v, then n=7) . The remaining differences 
"between the observed and the calculated values are distributed 
as will be later shown, almost equally to the members applied 
for m greater than 4 and that for n greater than 6, so that the 
two limits may be considered equal in a practical respect. 

Since the calculation cannot be presented completely, as [23 

already remarked, I dispense with it in the hope that I can provide 
the complete information later, even in an abstract, and now turn 
to the results. 

Prom the normal equations, formed in the previously indicated 
manner, the coefficients of the series for aX sin v, 3 Y sin v 
and YZ result (under the common designation pjj, qjj for reasons 
of brevity), compiled in the two columns designated with IV 
and in the last column of Table VII (the significance of the 
numbers in rows 1, II, III will be explained on page 38 ). These 
same numbers are found again in Table VIII, rounded off to whole 
numbers of the unit employed here and in somewhat altered, easily 
understandable arrangement. The fact that some of the numbers 
occurring here still contain decimal places is a result of the 
conditional equations (6). The values compiled in VIII are assumed 
as final and employed as a basis for all further calculations. 

It is not necessary to explain that the small deviation 
introduced through the rounding off causes only a negligible 
increase in the error square sum as a result of compensation of 
the best values so that rounding off is thoroughly permissible; 
in contrast, it could be considered unnecessary, since further 
calculations have no noticeable advantage because of this. Prom 
this standpoint, the unaltered values given in Table VII would 
be justifiably preferred. The alteration chosen by me, however, 
has another purpose. The numbers found in VIII without decimal 
places are not to be considered approximation values, only rounded 
off one place more than those in VII, but rather they represent 
precise values, exactly defining a definite possible state of 
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the magnetic field of the earth, although only approximating 
the real condition. The numbers provided with two decimal places 
do not contribute to the definition of this state, as they result 
from the other coefficients on the basis of the conditional equ- 
ations (6). The values indicated here are therefore to be consid- 
ered as rounded off, and can be calculated more precisely with 
the aid of those equations (in the numerical calculation of X 
from the series for ax sin v, it is not practical near the 
pole, and impossible at the pole itself to form the sum 


a sin v 


(B» J*S + + BlBl . . .) odor BJ JjjL + _*L 


4- 


since this is converted into an indefinite expression for v-0° 
and for v=l80°. The conditions are the same in the calculation 
of Y. In order to avoid this negative situation, the series /26 
is converted, taking into consideration the conditional equations 
existing between the coefficients (expressing the existence of a 
clearly defined horizontal force at each pole) , possibly carried 
out in different manners. It is easiest to select the form 


ST* 


asmt> 






-pirn 


ft 


a suit) 1 ” a sin t) 

05" : a ° 0 — erg", 


Bl 




a sin v 




:C *-+ J 1/3 


because then the same function tables may be employed, no matter 
how far the derivation is carried. In consideration of this, 

I employed the first two coefficients (with n=m and n=m+l) as 
functions of the following in the case of m=0 and then to retain 
the same form, also in the remaining values of m. It might be 
preferred from another point of view to proceed in reverse and 
express the last coefficients of each series by the previous ones; 
however, the procedure is a question of subordinate significance.) 

The same considerations, leading to the determination of weights 
of the equations (p. 22 , 23 ), make the derivation of medium or 
probable errors of the calculated coefficients (in VIII) from 
the differences present between observation and calculation appear 
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insignificant as long as these differences may he substantially 
reduced by a somewhat greater expansion of the series derivation. 

I prefer to provide the differences individually in order to permit 
clear judgement on the usefulness of the analytical representation 
gained, but then to undertake a mere estimation of the average 
error of the coefficients on the basis of the average uncertainty 
to be set for the observed values. Beth data together show whether 
and which expansions of the series development are necessary 
in the present state of our empirical knowledge, when precision 
of the representation adapted to this is to be attained, and which 
further improvements may only be gained by new certain observations. 

The following values of the coefficients k, K 1, L m, M 
from 5 to 5 degrees for the entire range of v~0° ’to v= 180° in the 
tables XI a, b, c are found from the numbers of table VIII. 

The numbers are indicated up to a tenth of the unit applied here 
so that later new calculations may be supported on precise dif- 
ferential values. The deviations of these calculated values [2 

from those listed in VI a, b, c are found in the following tables 
XII a, b, c, rounded off to whole numbers so that they are more 
understandable. All these differences and those following are 
formed in such a manner that they represent the excess of the 
observed value over the calculated value. They therefore designate 
the portion of the observation results not yet reproduced by the 
analytical expression and are of the same type as these. Any 
improvement in the latter may not be directly considered as cor- 
rection in the differences, and when these are represented completely 
or partially analytically, the results are to be added without 
alteration to the already determined expressions. 

The values designated as observed for the coefficients k, K 
etc. are only indirect, in contrast to those derived from the 
spherical function series. Those numbers can be termed actually 
observed, defining the empirical facts employed as a basis, i.e. 
the values of X , Y, Z in the 1800 points of the earth’s surface 
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described in the beginning (the fact that these are also already 
more or less the result of processing the direct measurement results 
is not important for the present report) . A conclusive judgement 
can then be founded only on the comparison of the values determined 
at these 1800 points according to Dr. Neumayer's charts with those 
following from the representation in Table VIII. The tables 
in XIV a, b, c serve for this comparison, but in these only the 
1 56 points already occuring in Table V — from 10°to 10°in 
latitude, from 30 °to 30 °in longitude — have been taken into 
consideration, sufficient for gaining a survey. Each of these 
difference tables is divided into three sections. The first, 
headed by A ' , contains those differences, occurring because of 
the limitation of the trigonometrical series derivations for the 
members of the first four orders, in other words, the excess 
of the observed values (in V) above those calculated according 
to the tables VI a, b, c. The differences produced by the limi- 
tation of the development according to spherical functions are 
found under the designation A ' * , resulting directly from the 
numbers in Tables XII a, b, c. The sum of those differences, 
i.e. the total difference between observation and calculation 
finally forms the third section A , of the fables. 

In order to evaluate the numbers given here, it must be noted 
that the average intensity of the total force of the earth's 
magnetism is approximately equal to 5 0000 (with the extreme values 
of 26000 and 70000 ) , measured in the same unit, while the average 
values of X, Y, Z with the exception of the signs are applied at 
approximately 3°000, 6000, 4000. A direct comparison with the 
deviations obtained in earlier calculations of potential is not 
possible, because these can only be expanded to the spherical 
functions of the fourth order (at least insofar as they have been 
published). A comparison may be gained indirectly, however, 
with the aid of my provisional calculation mentioned in the begin- 
ning, also proceeding to the fourth order. In this case, the dif- 
ferences A were almost 1.8 times as large as the definitive calcul- 
ation presented here, although only the component An was reduced 
in these. Specifically, the average differences calculated app- 
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roximately with the aid of Fechner's equation from the averages 
of their 1 56 absolute values 


in the case of x T 

after the first calculation ± 762 610 

after the second* calculation ± 336 434 


z 

1331 

819. 


The introduction of the spherical functions of the fifth and 
sixth order (j.n the case of aX sin V those of the sixth and seventh 
orders) therefore produces a substantial improvement of the results 
Of the calculations carried out earlier, the Erman-Petersen calcul- 
ation is best suited for comparisons, because only in this case 
are the differences between observation and calculation given with 
respect to the components of force These differences, however, 
refer to other points than those chosen here and their number 
is also less (90 versus 156) » however, this condition should 
have only a slight effect on the comparison of both sides of the 
results. It is more important that the 9 coefficients of the trig- 
onometrical series employed by Erman and Petersen are only derived 
from the elements of the 9 points of the corresponding parallel 
circle subsequently taken into consideration in the check, so 
that the differences designated with A* disappear there and the 
entire differences appear relatively too small in the case of 
Erman-Petersen. According to the data on p. 24 of the above- 
mentioned publication, the probable error of a calculated value 
now amounts in conventional units to 


in the case of-^ T z 

± 21,41 15,98 29,10, 

producing the average error in the unit employed here through 

n 

A. Erman and H. Petersen, "Fundamentals of the Gaussian theory 
and the phenomena of the earth's magnetism in 1829, " Berlin 
1874, p. 24, 


multiplication with 24.941 : 0.6745 
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in the case of X y z 

db 1109 828 1607 


These differences are 1*3 times as large on the average 
as those given previously, corresponding to my provisorial calcula- 
tion and 1.9 times as large as the actual differences A" to be 
compared with them. It cannot be determined with certainty even 
in an intensive study how far this result, favorable for the 
latter calculation, is based on the fact that the new observation 
material is better than the old material. Partially, the superiority 
of the new calculation is doubtlessly based on the fact that each 
component is represented independently in it (of course, taking 
into consideration the purely analytical, non-physical conditional 
equations combining X and Y), while all three components were found 
according to the old method, on the basis of certain physical 
hypotheses (p. 3) by a mutual compensation. Therefore, the number 
of the constants to be determined and employed as a basis for the 
representation is three times as great as here. Therefore, 
the new calculation must already produce slighter differences, 
when these are based only on observational errors; this must be 
the case to a higher degree, when the theoretical assumptions 
of the earlier calculation are not fulfilled (the hypothesis 
of the existance of a potential of exclusively internal forces) . 

The resulting contradictions between the values of the three com- 
ponents not fitting into the assumed form must be included in the 
final differences, so that these cannot be reduced below a certain 
amount by any expansion of the series, no matter how far. In con- 
trast, the differences in the independent calculation of individual 
components, as appearing in XIV a, b, c may be made randomly small. 
The superiority of the new method over the old must therefore 
become more emphasized as more members are taken into consideration 
in the derivation. Previous experience corresponds to this ex- 
pectation and therefore supplies an indirect proof for the fact 


that the physical prerequisites * of the old theory could not be 
strictly fulfilled. While the addition of the spherical functions 
of the fifth and sixth order in present calculations resulted in 
a substantial reduction in differences, as was already noted (p. 28), 
the same expansion of the derivation in Dr. Neumayer’s calculation £ 30 
of potential (the only one expanded up to now beyond the fourth order) 
had no substantial use 1 ), although in this case, the members dependent 
on 51 and 6X were taken into consideration, reducing the differences 
A* in contrast to the present calculations. 

Although the result of the new calculation may be designated 
as relatively favorable, it is still not satisfactory. The diff- 
erences compiled in the tables XIV a, b, c are apparently still 
substantially greater than the possible errors of the observed val- 
ues (i.e. taken from the charts), for which it is hardly possible 
to gain even a relatively certain estimation. Therefore, a contin- 
uation of the development could be considered. Some considerations 
also at least contradict a considerable expansion of the series 
and make it appear useful to proceed in steps, empirically deter- 
mining the most practical limit for the derivation on the basis of 
successive results. The larger the number of considered series 
members chosen and therefore the unknowns to be considered, the 
greater is generally the influence of observation errors on the 
coefficients. A continuation of the derivation therefore causes 
better adaptation to the finite number of observation values avail- 
able, but the entire representation becomes less certain, worsening 
in the intermediate points. This is to be anticipated in especially 
high degree in the case of the calculation basis employed here, 
containing two large areas completely free of observations, the 
two polar circles of 30°spherical radius. When the precision 
is too great in the representation of the values given for the zone 
from 60° to the north to 60° to the southern latitude, the derivation 
for the poj.ar areas would be completely unuse able., In addition, 

) 1 

Neumayer, "On magnetic surveying of the earth," lecture given 
at the meeting of German Scientists and doctors in Bremen, 1890. 


S*) 


there is the consideration that an almost complete representation 
of the observed values through an analytical expression is not even 
the most important task. Especially because this task can be solved 
theoretically with random approximation, as it is expressed here, 
while the lack of rules for the individual phenomena excludes rapid 
convergence of the series or a regulated correspondence of coef- 
ficients, the actual performance of the series derivation up to 
members of very high order is theoretically uninteresting and prac- /31 
tically worthless. The main significance of the analytical represen- 
tation is that it alone can supply a reliable basis for the 
physical theory of the phenomenon, that even the first question 
about the spatial origin of the force must remain unanswered, f-jv 
example, without it. The analytical representation permits individual 
consideration of the particular solutions of the fundamental dif- 
ferential equations of the problem by separating the entire picture 
of the phenonemon without rules into a (infinite) sum of regular 
components. The result is that the conclusions to be drawn from 
the series derivation are dependent on the expansion of the series. 
Therefore, the number is not as important as the most precise pos- 
sible determination of the series members to be calculated and the 
most important task must be presently considered an essentially 
precise determination of the first coefficients of the series. 

Of course, where each normal equation contains several unknown co- 
efficients as a result of the distribution of the observed values 
so that these cannot be obtained independently of one another, 
both cannot be completely separated. Up to a certain unpredictable 
limit, therefore, only the expansion of the series derivation 
leads to a more precise determination of the individual coefficients, 
especially the first and most important. A completely independent 
calculation of each coefficient and the possibility of obtaining its 
value free of any uncertainty attributed solely to observation 
errors requires the knowledge of the functions to be represented 
(for the components of force) on the entire surface of the earth, 
therefore appearing desirable here once again. 
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In the proceeding observations, only the purely analytical 
representation of the empirically provided distribution of force, 
clearly defined by the numbers of Table VIII have been treated. 

With the aid of the equations compiled in the first sections, 
a second equivalent representation is to be derived, determined 
in form according to physical considerations, serving as a basis 
for studying the causes of this distribution of force. Por this 
purpose, the functions U and W are to be derived first from aX sin v 
and 3Y sin v. Their coefficients are found in Table IX, also con- 
taining those of (Vi b) , i.e. the arithemetic mean of U Q and W Q . 

The circumstance that U and W are not in agreement expresses accor- 
ding to what has been said earlier that the horizontal forces 
in the earth's surface cannot be attributed completely to a poten- /3 2 
tial. When an electrical current is assumed as cause for this, 
penetrating the earth's surface perpendicularly* the intensity 
in relation to the surface unit can be calculated, in other words 
the current density i with the aid of the equation (14) from the 
difference (lnUU)* An enclosed derivation results for the almost 
proportional quantity a3bi, with coefficients to be found in Table 
X* The average amount of bi on the entire earth's surface results 
from this at about 110 t i.e. 0.0011 cirf u g w s~^. Since b, the 

Q 

polar radius of the earth, amounts to 6.356x10 cm, the average 

—12 A „.•] 

result for i results with 1.7x10 cm * g a s or approximately 
1.7x10 Ampere i cm . (Of course, square agents are meant here, 
since the algebraic sum of the currents and the simple average of 
the intensity values disappear.) Therefore, on an area of 1 square 
kilometer there is an average current of 1/6 Ampere. The physical 
possibility of such a weak current cannot be argued j there is no 
reason from this side to doubt the actual difference of W and U. 

Another question to be 3tudied is whether or not the uncertainty 
of the coefficients of U and W is so large that the difference 
of the two functions becomes illusory. 

Prom the potential V of the main portion of the horizontal 
forces, chosen as large as possible, and from the vertical component 
V, the potentials V^andVgOf the inner and outer agents result 
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according to the equations 15/18* The coefficients of both, 
or rather the proportional quantities (V^i b) and (V &| b) are given 
in table X. It can be seen that v_ isactually very small, as was 

C* 

predicted according to the earlier potential calculations; the 
average amounts to only approximately l/40 of that of V^. Of course 
even here, the certainty of results must still be examined. When 
I first dispense with this, I can express the result according 
to the calculations in the following manner. 

The empirically determined distribution of the magnetic forces 
of the earth in the earth's surface, as they were represented 
in Dr. Neumayer's charts for the moment 1885,0, is based mainly 
on causes situated in the earth's core, but cannot be attributed 
exclusively to these. A small portion of the force, (about l/40 
of the entire amount) can be attributed to causes outside of the 
earth's surface; a further still somewhat greater portion (especially 
the horizontal component) points to electrical currents penetrating 
this surface. 


The two representations of the magnetic state of the earth, /53 
defined on the one hand by the numbers in VIII and those in X on 
the other hand, have this same value, as was already noted. One 
can be derived from the other, taking into consideration the earlier- 
mentioned side conditions (p. 10). It is therefore understood, 
but still deserves special emphasis, that the representation given 
in X deviates from the observed values only by the differences 
already present earlier (in the case of VIII), expressed in XII 
and XIV. It is therefore possible in principle to have the rep- 
resentation through V^, V a and i approach the actual, empirically 
determined state with random accuracy. 

The currents parallel to the earth's surface can easily be 
calculated from the values found for V^ an&V a (compare equations 
22) , producing the forces expressed therein. I will dispense with 
providing the coefficients of the corresponding derivations, since 
they may be considered relatively uncertain. In the case of V^, 


I 
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the possibility that at least a noticeable portion is to be attri- 
buted to a direct magnetization of the earth's crust cannot be 
rejected; on the other hand, in the case of V a the uncertainty of 
the values found, although not absolute, but still relatively large, 
is too considerable to provide sufficient value for the derivation 
of detailed conclusions. Therefore, I will limit myself to a 
summary, only characterizing the magnitude of the currents necessary 
for explaining the phenomena. In this case, I imagine in the 
manner described earlier (p. 14) that the entire flow is compressed 
through vertical shifting into two infinite, thin layers infinitely 
adjacent to the earth's surface on both sides and indicate the 
current, passing through a cross-section perpendicular to its dir- 
ection of 1 cm in width. The average intensity expressed in this 
manner is rather precisely l/2 AMP: cm for the internal and 

approximately 1/700 AMP: cm for the external current, when the 

meridianal component, rather slight in the case of V. and hardly 
certain in the case of V a , is disregarded. The direction of both 
goes from east to west. 

It is now appropriate to list some further results of a similar, 
general significance, the moment and the axial direction of the 
earth magnet, determining the effect at a distance. These two con- 
cepts must be limited to the portion of force originating within 
the earth. The external agents expressed in V a could.be represented 
by a moment of defined direction and size only with specific assump- 
tions on distribution, when this concept is understood in the custo- /54 
mary manner. In contrast, however, a vector corresponding to the 
moment could be indicated in the case of VI in another interpretation, 
determining direction and strength of the magnetic force at the 
center of the earth. This data would probably not be especially 
interesting; therefore, I will not provide the appropriate results. 

The value results for the magnetic moment ( it could be termed 
that of the solid body of the earth) under the given limitation 
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M = 8,3481 • 10" cm* g\ s~ J (log M = 25,92159) 


and as axial direction those according to the point 


t> = 168°32',1 X = 111 0 29', 4 

or u — 168° 34', 3 X = 111 0 29', 4. 

The axis therefore connects the point of 78°34',3 northern 
latitude and 68°30',6 western longitude from Greenwich with that 
situated at the 78°34'3 southern geographical latitude and 111°39',4 
eastern longitude. 

For comparison, I mention that the moment according to the 
calculation of^Neumayer and Petersen amounts to 0.32237R^, i.e. 
8.3324xl0 2 ^ cm* 3 g’ a ’ s-^C with R= 6 . 370x10® cm) and that the axis 
is directed to the point (f>=-78 0 20', X =112°43 , » The slight dev- 
iation of the numbers from those given previously shows clearly 
how much the causes of the magnetic force of the earth is situated 
in the earth's core. The question still remains to be studied 
on the certainty of all these results and which conditions are to 
be fulfilled to increase the certainty. 

When the indicated values, considered as observed, were com- 
pletely free of error* i.e. when they correspond in complete precision 
to the actual space of the earth's magnetism in a defined moment 
(even with elimination of exactly defined local and temporal dis- 
turbances) , the differences (portions of which are given in Tables 
XIV) of the calculated and observed values would be a consequence 
of the early interruption of the series and it would be possible to 
reduce them through a continuation of the derivation. It would 
not be practical, however, even in this favorable case, to continue 
the derr at ion so far that the differences disappear; this would 
only be advantageous if it were known a priori that all members 
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of higher order in the exact representation (only obtained by inte- 
gration over the entire surface of the earth) , not determined 
in the calculation of individual values, are actually negligible, / 55 
If this is not known — and in the present case the contrary can 
be concluded from a mere survey of the entire irregular and slowly 
decreasing coefficients — the representation of the above-mentioned 
values would be improved, as was noted earlier, but those of the 
entire state would be negatively affected. This is the case to 
an even greater degree, when the given numbers are erroneous, so 
that the task assumes the character of a compensation calculation. 

Therefore, as long as the observation basis cannot be expanded, 
a type of scanning procedure is recommended, providing an evaluation 
through successive expansion of the series derivations on which 
alterations are undertaken for the first, most important coeffic- 
ients and the degree to which their terminal values may differ from 
the limit they approach in the further derivation. Now I have 
undertaken a provisorial representation, proceeding only to the fourth 
order, as was noted earlier? moreover, I have also carried out 
a greatly abbreviated derivation expandeu only to the second order 
(in the case of aX sin v of course to the third order). (Since 
the members of even and uneven order are calculated independently 
of one another, so that the two intermediate cases are processed 
for which the derivation goes to the third and the fifth order.) 

I shall not provide the results with the exception of some examples, 
but only give a summarizing result that those portions of the de- 
rivation dependent on a multiple greater than double the geographical 
length X may only be assumed as raw first approximations, and, on 
the other hand, that the first coefficients are to be determined 
so precisely for m=0 and m-1, to a slighter degree also for m=2, 
that a continuation would not alter them very much. Careful ex- 
amination of the differences given in Table XII leads to the same 
conclusion. A further circumstance is important for the evaluation 
of results. The two derivations apply to aX sin v and 3Y sin v 
are connected by a number of purely analytical conditional equations 
not proceeding for physical reasons, according to earlier explanations 
(p. 11). It is now shown that these equations generally introduce 
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a surprisingly strong constraint. I first determined the coeff- 
icients of both derivations independently and subsequently calculated 
the alterations caused by the introduction of the conditional equa- 
tions and added these « The result of these alterations was very 
considerable in almost all cases, as was already mentioned; they 
decrease with the continuing expansion of the series derivation 
and are smallest in the first members of the series. The main 
reason for the contradiction appearing in these differences lies 
in the irregular distribution of the magnetic force of the earth, 
only represented by slowly converging rows; the lack of observed 
values from the polar regions has a distorting effect especially 
here, and this is easily understandable when it is noted that those 
conditional equations refer to the arrangement of forces at the poles. 
The errors involved in the given values have only an indirect ef- 
fect, insofar as they supply relatively large amounts in their 
irregular distribution, especially to the members of higher order 
numbers. 

An example is the portion multiplied by cosA of the series 
for aX sin v and the related factor of sinA in the series for 
3*Y sin v. 

3 

inlhen the derivation is expanded up to in. the case of X, 
that is R£.in the case of Y, the coefficients defined without consider- 
ation of the conditional equations are 


— 1943,9 198,9 732 , 9 . 

“ 1341,3 1273 , 8 ; ' 


while the values actually to be chosen and satisfying the conditional 
equations are as follows t 



rite 1 


— 1872*1 — 455,5 1128,5 

- 1018,6 1051,2. 
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When the derivation is carried two steps further, the numbers 
result as independently calculated: 


— 1930,5 205,7 757,6 — 997,5 178,9 

— 1290,8 1272,5 — 525,9 189,7 


and when compensated: 


— 1913,2 270,8 817,1 — 790,2 308,7 

— 1312,3 1249,6 — 644,1 116,8. 


When two further members are taken into consideration in each 
series, the corresponding coefficients are 


— 1932,3 

— 1279,8 


270.4 752,1 — 899,9 172,0 

1272,1 — 505,5 189,1 147,9 


313,4 —18,9 

~ 1,7 


and 


— i920,9 258,4 " 788,5 — 924,4 226,9 . 273,9 5.4,6 ;. 

1276,9 1260,9 — 495,6 166,1 169,7, — 38,7. , 
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Consideration of these numbers demonstrates that only the first 
coefficients are somewhat certain and that a separate calculation 
of the individual coefficients is desirable. Such a procedure, 
however, requires the knowledge of the magnetic elements in the 
polar areas (although they are not undertaken through evaluation 
of surface integrals but according to the Neumann method) . 

When the corresponding parts of the series for (V : b) and 
a$bi are calculated from the previously reported partial develop- 
ments of X and Y, similar differences must of course result in 

■) 

the progressive expansion of the series. The coefficients of 

cosA , R^cosA .... in (Vj b) in the three cases under consideration 

are 

— 1019,7 1116,3 ' * ’ ' l “ . • 

— 1506,6 1252,4 — 591,8 137,0* / 

— 1220,4 1231,6 — 451,3 156,9 140,4 —9,8 


and those of R^sinA , R^sinA ... in agbi 


— 45,8 • . . 

— 30,4 — 52,1 — 26,0 •* \ * • 

— 33,3 3,3 — 39,7 47,9 — 54,5. - 


The coefficients of the series for yZ limited by no conditional 

equation (with the exception of the one Jq- 0) demonstrate generally 

small alterations in the continued derivation. The coefficients 
1 2 

belonging to R-^cosA , R^cosA ... will be listed here as illustrations. 
In the three cases they are expressed as follows: 

3083 — 3819 * \ * V - , 

. 2893 —3813. 1984 — 868 

2847 — 3841 1900 — 909 —607 — 132.' 
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Concerning the above-mentioned conditional equation 
this is fulfilled almost without any constraint, in contrast to 
the previously mentioned equations « When the coefficients 
jg, jg, Oq are determined without consideration of that equation 
the results are 5*9 — 708. 5 — 1265*2 — 32 * 7 , while the 

results under consideration of the conditions are 0 — 7 ° 1 *^ 

1272.6 — 40.0, The previous considerations give the impression 

that the relatively small differences for concluding the existance 
of a current penetrating the earth's surface and of a potential /38 
of external agents must be termed illusion. This impression is 
intensified by further considerations » Table VII contains four 
columns each in the case of X and X, but only the last one has 
been mentioned up to now. The numbers in the three preceeding col- 
umns designate the coefficients of the series under the following 
assumptions s 

I. i«0, V =0, i.e. the magnetic force of the earth is based 

El 

exclusively on the effect of inner causes, as was assumed by 
the Gaussian theory. Under this assumption the series for 

aX sin v and pY sin v are calculated from those for YZ. 

II. The series are derived without consideration of the 
conditional equations (some of the numbers listed here 
were already mentioned (p. 36 )). 

III. i=0, i.e. the series has been calculated under the 
assumption that the entire force exhibits a potential 
partially based on external causes. 

It can be seen that the result of the calculation carried out 
without any hypothesis, listed in column IV, does not follow any 
better from II than does III, while I deviates somewhat more. 

In this case, it must be considered that I is supported on the values 
of yZ , assumed to be unaltered; in a common compensation, of 
X, Y, Z (as is undertaken in the Gaussian procedure), instead of I 
a value system situated between I and III was the goal. 
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In order to illuminate the matter from still another side, 

I have calculated the values of ^K-jl-jI^, although only for the por- 
tion dependent on cosX and sinX and have formed their differences 
to the observed values. These are found in Tables XIII a and b. 

It should not be forgotten in comparison that the error square sum 
is not formed from these values themselves, but from their product 
with sin u (or strictly with a sin v in the case of X, $ sin v 
in the case of Y), so that larger differences in higher latitudes 
have less effect than such near the equator. In fact, it can now 
be seen that the differences in the case of IV more closely approximate 
those in the case of II than the other two, of course representing 
the error minimum. This supports the fact that the representation 
of the actual distribution of force exhibits clear progress in 
spite of all uncertainty of individual components given by their 
totality, when the hypotheses that i and V 0 disappear are dispensed 
with. 


I have not yet taken into consideration that uncertainty /59 

stemming from the errors in observed values. The preceeding con- 
siderations are hardly affected by this, since the effect on the 
four parallel calculations cannot differ very much. The amount and 
direction of this effect on the final results will now be examined. 

Of course, it is presently impossible to carry out precise calcula- 
tions and these will not be possible for a long time, because 
they assume the knowledge of extremely alternating uncertainty 
for the individual evaluated observations and surveys as well as 
those attached to the assumed secular variations, and while it 
is still difficult for broad areas to gain any useable values 
of the magnetic elements, it is much more difficult to form a 
judgement on their approximate precision. Even a raw estimation 
is not possible in this case in a direct manner, since it would 
have to return to the original material processed in the charts. 

A conclusion about the possible amount of average errors of the 
observed values can be drawn indirectly, however, from a consideration 
of the differences between the observed and calculated values. 
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The slow reduction in coefficients of the series development with 
increasing order number n removes any doubt that the remaining 
differences are attributed to a great part to calculation of the 
insufficient expansion of the series and the grouping of these dif- 
ferences still far removed from irregular arrangement leads to the 
same conclusion. When it is now observed how substantially the 
average deviation was reduced according to earlier data (p. 28) 
by the addition of the members of the fifth and sixth order, it 
can be assumed without a doubt that the average amount of deviations 
attributable to the inaccuracy of the observed values is at most 
half as great as the entire amount of these. Therefore, I will 
assume in round numbers that the values of the components of force 
have an average error of i 200 in the case of X and Y and of 
400 in the case of Z for each of the 1800 points utilized and these 
numbers can probably be considered an upper limit. Average errors 
of -20' and ^0.002 cirT^" s“ 1 could be considered approximately 
equal in value in the declination and horizontal intensity and er- 
rors on the order of magnitude of 1° in inclination. Of course, 
all unknown effects of temporal or local disturbances, the inaccuracy 
of the observation itself and the uncertainty of reduction to the 
normal epoch, finally also those in the drawing of the isomagnetic 
lines and in the removal of values for certain points between the 
poles from the charts are contained in these numbers. 

Since the extreme differences existing between the individual 
areas of the earth's surface with respect to the reliability of our 
knowledge on the magnetic state is completely neglected in this 
procedure, it can only be designated a summarizing estimation; 
however, substantial distortion of the conclusions to be drawn from 
this must not be anticipated, since the only consideration is the 
size of expected average uncertainty. * 

The coefficients k m ,K m . . . of the trigonometrical series are 
derived from 36 individual values each; therefore, the average 
error results from the above-applied amounts by adding the factor 
1/6 or l/6 according to whether the index m is equal to 

0 or not. When e is the absolute amount of error, then 
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C V 

a fflr m = 0 b bei X und T ; f = 83 D bei Z: e = 67 

ftlr m = 1, 2, 8, 4 bei X und Y. : e = 47 bei Z : f = 93 


Key: a. For b. In the Case of c. And 

( A justifiable further rounding off of these numbers will not be 
undertaken in order to prevent distortion of the relationships 
of various values.). The normal equations for determining the 
coefficients B, C, D, E j , k from the values of k m) ... also lead 
to the calculation of the average error of the former in a known 
manner. I again provide only some results , sufficient for character 
ization. In this case, I would like to note that I have assumed 
the values given for e as average error of the quantities ctk sin v, 
aK sin v. . . in order to substantially shorten the calculation 
for X and Y, although they belong to k, K... The errors of B, C, 

D, E resulting in this manner, a few of which will now be listed, 
are therefore somewhat too large, but this is not a problem since 
the only consideration is the calculation of an upper limit for this 
error. 

In the case of m=0, the following anticipated average errors 

of the coefficients B n , D n , J n result 

o o o 

n = 0.1 2 3 4 5 6 7 

- a bei X ± (90) (44) 15 13 19 13 17 14 

bei Y ±(90) (24) 15 11 19 3 17 . 

bei Z ± 52 51 60 45 55 ,36., 

Key: a. In the Case of 
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The values related to B“, B?, D“, Clare enclosed in parentheses, 
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while these coefficients are completely determined by the others 
(p, 23) and therefore have no independent effect on the values 
of the function presented. In order to characterize the effect 
of the greater or smaller expansion of the series development 
on the average error of the coefficients , insofar as these stem 
from errors in observation material, a typical example may be suf- 
ficient. According to whether the calculation of the coefficients 
0O» 0o» Oo is interrupted with the first, the second or, as was 
undertaken in the present report, with the third, the average 
errors are 


Us± 


0,24 <; 


0,37 e, 0,37 <; 


0,77 e, 0,67 t , 0,54 e. 


These numbers demonstrate that an uncritical continuation 
of the derivation without expansion of the the empirical basis 
would not lead to an improvement, but a worsening of the results. 

For m=1 , which I will describe in detail later, since the re- 
sults for m s 2, 3 1 4 are on the same order of magnitude, the following 
amounts of error result 



n 

•sz 

1 

2 

3 

4 

5 

6 

7 

bei 

X 

± 

8 

8 

14 

10 

19 

13 

21 

bei 

Y 


8 

8 

10 

11 

12 

14 


bei 

Z 

± 13 

14 

13 

17 

16 

21. 



Key: a. In The Case of 

These error amounts depend of course on the formation of the 
normal equations; it would be decreased still more if the obser- 
vations were distributed evenly over the entire surface of the 
earth . 
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The average errors given in the preceding, moreover representing 
upper limit values, are slight and hardly need be considered es- 
pecially in the first coefficients of the series. The circumstance 
is a result of the fact that the derivation is based on the elements 
of an extremely large number of points. The analytical represen- 
tation of the components of force can therefore be designated 
as very precise and reliable with the exception of the systematic 
error caused by series interruption. 


The series for V\ , V & and i are much less favorable. The an- 
ticipated average errors are not only relatively larger in this 
case, but they achieve the maximum precisely in the first and most /42 
important coefficients. I will limit myself here to the presentation 
of an example? specifically, I chose those coefficients of the 
above-mentioned functions for this purpose, resulting from the 

previously indicated coefficients for X, Y, Z. These are the factors 
12 12 

connected to Rq, Rq**' and those with R-j , R| • . . • I will compile 
them in an arrangement similar to that of the previous examples: 



m 

= 0 

n — 

1 

2 

3 

4 

5 

6 

in. 

F. 

bei V{\ b 

± 

31 

15 

12 

7 

6 

3 

m. 

F. 

bei V n : l 

± 

41 

16 

14 

r 

7 

3 

ra. 

F. 

bei a [3 bi 

± 

14 

7 

18 

7 

16 



Ml 

as 1 

n — 

1. 

2 

3 

4 

5 

6 

m. 

F. 

bei V t \ b 

± 

53 

27 

33 

19 

23 

16 

m. 

F. 

bei V.ib 

± 

75 

30 

35 

20 

23 

16 

m. 

F. 

bei a [ibi 

± 

38 

15 

22 

11 

15 



Key: a. Average Error in the Case of 

At least in the case of m=l (and it is similar in the case of 
m=2, 3, 4) these error amounts, however upper limit values, often 
exceed the values of the coefficients of V a and i. Especially these 
functions, .interesting in a physical respect, therefore become 
so uncertain that it appears doubt ful whether they are perhaps 
only a product of observational errors,, Some early experience, 
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however, supports the fact that this is not the case (p. 29 » 38) 
and the individual derivation coefficients may indeed he very un- 
certain while the function presented assumes values of much less 
uncertainty at some positions. This is caused hy the fact that 
the coefficients are dependent on one another, since they partially 
depend on the same variables (coefficients of the series for X, 

Y, Z). 


The reality of V & and l are supported by the circumstance 
that at least a portion of these functions noticeably exceeds 
the error limits, specifically those belonging to the index m=0 
and therefore independent of the geographical longitude, so that 
it may not be considered a result of mere internal contradictions 
of the data provided. A more precise consideration intensifies 
this impression. The component free of x in a(8 bi follows from the 
portions of the development of (3Y sin v also only dependent on 
v with coefficients designated on the individual parallel circles 
by 1 q. The condition that in the derivation according to X complete 
circles with evenly distributed observation values are present 
(in contrast to the derivations according to v, supported here /43 
only on the arc from v=30° to v=150°) has the effect that the coef- 
ficients of the trigonometrical series may deviate from the true 
values almost solely as a result of the errors in empirical data. 

The effect of these random errors on 1 Q was recently found (p.10) 
at an estimated ±33- When this is now compared with the values 
of 1 q given in Table Ylb; growing even in the relatively well studied 
medium latitudes of the northern hemo sphere up to 150 » on the south- 
ern hemosphere even up to the double amount, it can be seen further- 
more that these numbers exhibit a very clearly marked course, 
so it appears thoroughly impossible to attribute these to errors 
in observational material. This would only be possible when 
large systematic errors could be distributed over large areas 
(entire continents or oceans), a rather improbable assumption. 

When the existence of vertical currents i has been proven 
with a great probability by means of the considerable values of 
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1 0 , there is no reason to assume that the coefficients found in 
the derivation of i with the indices m~l, 2, 3t ^ are actually 
0 and that the values found for this factor only stem from obser- 
vational errors. These coefficients, however, have been determined 
with so much uncertainty according to the above discussion that 
an evaluation of the series found for i could only supply values 
of a rather doubtful nature. Therefore, X will not provide any 
information on the calculations I carried out in this direction. 

If there are currents connecting the space inside the earth's 

surface with that outside, a potential of external forces in the 

earth's surface is necessarily proven. Independently, at least 

*1 

the strong value of the coefficient of Rq, and probably that of 
R^» not attributed to faulty observation, also supports the actual 
existence of V_ . 

ct 

Let us survey the essential results once again in context. 

It can first be stated that the initial task, that of the purely 
analytical representation of the observed phenomena, has been solved 
in an almost sufficient manner in the expansion of the series 
applied here. The remaining differences demonstrate that the addi- 
tion of the members, depending on a multiple of 5 and at most 
also those depending on a multiple of 6 of the geographical lengths 
X and those of the seventh and eighth order in the spherical functions 
of the polar distance v would certainly be sufficient to provide / 44- 

a representation of the observed values with such precision as is 
desireable with respect to the degree of its inherent certainty. 

On the other hand, it is certain that the outermost limit is achie- 
ved in such an expansion of the series, and perhaps already ex- 
ceeded, beyond which the increasing random error of the calculated 
coefficients cancel the gain in more precise representation and re- 
verse it. This problem, however, would be somewhat reduced when 
the already very large number of 1800 points, from which the elements 
were utilized, is increased or when the derivation is supported 
only on the graphic representation of the elements; however, 
a considerable reduction in the resulting average error cannot be 
expected, because the largest portion of the earth's surface has 
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not been sufficiently surveyed so that the errors in adjacent 
points are considered mutually independent and may be treated as 
random errors. 

Accordingly, when the total representation can be designated 
essentially as satisfactory, only the first coefficient of the 
series can be considered reliable. As a result, the series for 
V,, V a and i, derived from these coefficients, expressing precisely 
the physically interesting results, are already less certain when 
considered absolutely, because each of their coefficients is based 
on several coefficients of the original series. Moreover, since 
V a and i are relatively small, the possible errors in this case are 
extremely large when considered relatively, so large that the 
actual existence of the physical conditions expressed by them is 
almost placed in doubt. And although both a potential of external 
agents and a component of the horizontal force not attributable 
to a potential may be considered proven, it must be admitted that 
nothing reliable can be determined on the precise shape of both 
with the present basis. Especially detailed knowledge of this shape, 
i.e. the distribution of the functional values of V a and i over the 
earth's surface is manditory, when the question now arising about 
the nature of these phenenomena is to be treated. For example, 
how important would it be, only to note one possibility, when 
an analogy with the configuration of water and land or a connection 
to climactic areas was demonstrated in the arrangement of the posi- 
tive and, negative values of i. After what has just been said, 
it needs no justification to dispense with speculations on such £/£ 
possible relationships, as well as any physical explanation of the 
phenomena, although several reasons could be found for this. I 
will only the note the conclusions, resulting, on the one hand, 
from the slow convergence of the series and, on the other hand, 
from the circumstance that the first members from which the magnetic 
moment originates dominates all others considerably; furthermore, 

I wish to mention at least the mutual effect of the inner and outer 
agents on one another. To treat this subject in more detail or 
some others is not only forbidden by consideration of the above- 
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described conditions? it would also extend 'beyond the framework 
of the present report, only aiming at the most exact possible ana- 
lytical representation of the observed facts. 

Eke utilization of material available from the areas north 
of 60 degrees in latitude and determination of coefficients in 
a modified manner already discussed (p. 21) would probably make 
the results appear somewhat more favorable. A substantial advance, 
however, with respect to the certainty of results can only be attained, 
as can be seen from the discussion above, when the calculation 
is based on complete knowledge about the distribution of the mag- 
netic forces of the earth in the earth's surface. Of course, this 
requires, above all, gaining new reliable material from the southern 
polar region, but also undertaking new observations in numerous, 
scattered, accessible areas with special consideration of determining 
the secular variation. Too much precision in detail and a very 
detailed knowledge of the distribution of force in smaller areas, 
although this may be important for other, specific tasks, may be 
termed less significant for present purposes, insofar as it cannot 
compensate for the lack of observations at other points. 

The demand results for future research from the results and 
considerations above that the consideration of the magnetic phenomena 
of the earth should be placed in the foreground as a unified pheno- 
menon to a much greater degree than has been done in the past. 

The numerous individual studies of special tasks provided for 
research from a multitude of aspects should of course not be under- 
estimated or neglected? however, although these may be important, 
they only obtain a true value through the classification in the 
relationship of the phenomena and this research should not be 
pushed in the background because of that. The most pressing task 
in the near future is to obtain a unified picture of the entire 
state of the magnetic system of the earth in its various relation- 
ships on the path of observation, first only in general, character- 
istics and then later more detailed, and to follow its constant 
alteration. A planned, organized procedure is required for this 
purpose? this becomes even more necessary with less means for 


research to solve these tasks. 
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